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We introduce a method to experimentally measure the monochromatic transmission matrix of a complex 
medium in optics. This method is based on a spatial phase modulator together with a full-field interferometric 
measurement on a camera. We determine the transmission matrix of a thick random scattering sample. We 
show that this matrix exhibits statistical properties in good agreement with random matrix theory and allows 
light focusing and imaging through the random medium. This method might give important insights into the 
mesoscopic properties of complex medium. 
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The propagation of waves in heterogeneous media, espe- 
cially in the multiple scattering regime, is a very fundamental 
problem of physics with numerous applications ranging from 
solid-state physics and optics, to acoustics and electromag- 
netism [ 1 ]. The behavior or such media is usually described 
by their average macroscopic properties such as the scattering 
and absorption mean free paths, and via spectral or temporal 
properties, such as spectral resonances or confinement times. 
Some of those quantities can be obtained through intensity 
measurements using for example speckle pattern correlations 
|2] or coherent back scattering experiments OH). 

A deeper approach for the study of complex media lies in 
the transmission matrix (TM) retrieval. This matrix is a sub- 
part of the usual scattering matrix as defined in Q for in- 
stance. Within this framework, the Green's function between 
an array of sources and an array of sensors is recorded in trans- 
mission. The knowledge of the TM brings more fundamental 
insight on the medium. One can for instance extract from the 
TM the single and the multiple scattering component 0, the 
backscattering cone, and the field-field correlations |7 , 8 ]. The 
distribution of the singular values of the TM should also be 
related to diffusion properties and might exhibit coherence ef- 
fect beyond the diffusive transport regime such as weak and 
strong localization effects. 

Furthermore, from an operative point of view, the knowl- 
edge of the TM of a complex medium offers new and exciting 
possibilities. For instance, using Time Reversal (TR), which 
is a broadband version of Phase Conjugation (PC), it has been 
shown that multiple scattering, far from limiting wave ma- 
nipulation through a random medium, can in fact greatly en- 
hance it. Such approaches, which are based on the reciprocity 
and the reversibility of the wave equations 0, have proved 
very useful in various areas ranging from focusing to imaging 
or telecommunication, in acoustic lfT0l - IT2ll . electromagnetic 
(l3][T4l or even seismology lfT5l . 

For acoustic (resp. electromagnetic) waves, transducers 
(resp. antennas) are natural local receivers and emitters in 
phase and amplitude, and oscillation frequencies are com- 
patible with electronic detection. Therefore, the TM can be 
straightforwardly measured. In contrast, reconstructing the 
transmission matrix of a complex medium is still an elusive 



problem in the optical domain. 

Nonetheless some recent experiments have demonstrated 
that it is possible to manipulate light at a mesoscopic level 
in a complex medium in order to focus light through lfl6l or 
in ifTTl a scattering medium, as well as couple efficiently to the 
open channels of a disordered sample 1 1 8 ] . These experiments 
were made possible by the emergence of Spatial Light Mod- 
ulators, which allow to control the phase of an optical field 
over a million pixels, and were based on algorithmic proce- 
dures which only give a limited information on the physical 
behavior of the medium. 

In the present paper, using a spatial phase modulator along- 
side a full-field interferometric method, we demonstrate that 
we are able to measure the monochromatic transmission ma- 
trix of a random multiple scattering medium. We first detail 
the experimental setup of the TM evaluation method. In a sec- 
ond part, we start by showing that the knowledge of the TM 
allows one to reproduce the results of fT6lL and make the link 
with TR experiments fTTl[T4ll . Then we demonstrate the po- 
tential of our approach with the detection of an unknown ob- 
ject after propagation through the scattering media. Finally, 
we study the statistical properties of the singular values of the 
TM, and show that it follows a quarter circle law typical of a 
random wavefield. 

We define the mesoscopic TM of an optical system for a 
given wavelength as the matrix K of the complex coefficients 
km n connecting the optical field (in amplitude and phase) in 
the m th output free mode to the one in the n th input free 
mode. Thus, the projection E™ 1 of the outgoing optic field 
on the m th free mode is given by 

n 

Where is the complex amplitude of the optical field in 
the n th incoming free mode. 

The sample under study is an opaque 80 ± 25 fim thick 
deposit of ZnO (Sigma- Aldrich 96479) on a standard micro- 
scope glass slide. The experimental setup consists of a diode- 
pumped solid-state single longitudinal mode laser source at 
532 nm (Laser Quantum Torus), a modulation part to illumi- 
nate the sample with a controlled wave front and a detection 
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part to measure the optical intensity transmitted by the scatter- 
ing medium. The laser beam is first expanded then spatially 
modulated by a liquid cristal Spatial Light Modulator (SLM). 
This modulator is a twisted nematic liquid cristal on silicon 
device (Holoeye LC-R 2500). Choosing a suitable combina- 
tion of incident and analysed polarizations, we achieve an al- 
most phase-only modulation fT9l of the collimated incident 
beam with a 2tt modulation in phase and a maximum resid- 
ual intensity modulation below 10%. The surface of the SLM 
is imaged on the pupil of a 20x objective with a numerical 
aperture (NA) of 0.5, thus a pixel of the SLM matches a wave 
vector at the entrance of the scattering medium. The beam 
is focused at one side of the sample and the output intensity 
speckle is imaged 0.3 mm from the surface of the sample by a 
40x objective (NA = 0.85) onto a 10-bit CCD camera (AVT 
Dolphin F-145B). The speckle is stationnary well over the 
measurement time (several minutes). 



intensity in the m th output mode is given by 




SLM 



FIG. 1: Schematic of the apparatus. The 532 nm laser is expanded 
and reflected off a SLM. Polarization optics select a phase only mod- 
ulation mode. The modulated beam is focused on the multiple- 
scattering sample and the output intensity speckle pattern is imaged 
by a CCD-camera. L, lens. P, polarizer. D, diaphragm. 



We choose the size of the input and output independent pix- 
els to have a perfect matching in size between a pixel and a 
mode, in particular a CCD pixel has the size of a speckle grain. 
Thus, in our setup, the input and output modes are the SLM 
and the CCD pixels respectively, and the TM corresponds to 
the system comprised of both the scattering sample and the 
optical system between the SLM and the CCD Camera. From 
now on, we fixed both the number of controlled segments on 
the SLM and the number of subdivisions measured on the ob- 
servation windows of the CCD to TV = 256. 

Measuring a TM in optics raises several difficulties coming 
from the impossibility of having direct access to the amplitude 
and phase of the optical field. Generally, one can have access 
to the complex optic field using interferences with a known 
wavefront and a full field "four phases method" l20l . For any 
input vector, if the global phase is shifted by a constant a, the 
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where s m is the complex amplitude of the optical field used 
as reference in the m th output mode. 

Thus, if we inject the n th input mode and we measure 
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For practical reasons, we choose the Hadamard basis as in- 
put basis over the canonical one. In the Hadamard basis, all 
elements are either +1 or -1 in ampltiude, which perfectly fits 
with the use of a phase-only SLM. Another advantage is that 
it also maximizes the intensity of the received wavefield and 
consequently improves the experimental signal to noise ratio 
(SNR). For all Hadamard basis vectors, the intensity is mea- 
sured on the canonical basis of the pixels on the CCD camera 
and an observed transmission matrix K Q i) S is acquired, which 
is related to the real one K by 



obs 
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where S re f is a diagonal matrix representing the whole 
static reference wavefront in amplitude and phase. Ideally, the 
reference wavefront should be a plane wave to directly have 
access to the K matrix. In this case, all s m are constant and 
Kobs is directly proportional to K. However this requires the 
addition of a reference arm to the setup, and requires inter- 
ferometric stability. To have the simplest experimental setup 
and a higher stability, we modulate only 65 % of the wave- 
front going into the scattering sample (this correponds to the 
square inside the pupil of the microscope objective as seen in 
Figure [TJ), the speckle coming from the 35 % static part being 
our reference. S re f is now unknown and no longer constant 
along its diagonal. Nevertheless, since S re f is stationnary 
over time, we can measure the response of all input vectors 
on the m th output pixel as long as the reference speckle is 
bright enough on the considered modes. We will quantify the 
effect of the reference speckle and show that neither does it 
impair our ability to focus or image using the TM, nor does it 
affect the statistical properties of the TM. 

A good way to confirm the physical relevance of the TM 
is to use it to focus light on any desired outgoing mode. In 
El H3|, it was demonstrated that using time reversal, one can 
take advantage of multiple scattering to focus on tight spots, 
thanks to the reversibility of the wave equation. In this con- 
figuration, the Green's functions between the array of sources 
and a point at the output of the disordered medium is sent re- 
versed in time, and the waves generated converge naturally 
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towards the targeted spot. TR being a matched filter ifTTTl . 
the energy is maximized both temporally and spatially at the 
intended location. The monochromatic equivalent of TR is 
phase conjugation, which can be straightfully done using our 
acquired TM. We expect similar focusing results as the ones of 
Vellekoop lfT6l . which were obtained through a procedure that 
ensures an optimum phase modulated wavefront maximizing 
the energy of a given mode. Denoting E tar9et the output tar- 
get vector, the input vector that approximates the desired pat- 
tern at the output for a perfect phase conjugated focusing is 
given by 

E in = K \ garget ^ 

where t stands for the conjugate transpose. Thus, the theo- 
retical effective output vector is E e ^ = o^ oc E target where 
Qfoc _ [ s ^ e so called time reversal operator |2T|| . The 
result of the theoretical focusing on the m th outgoing mode is 
the m th column of 0^ oc . Since our setup is limited to phase 
modulation only, and given the fact that we do not acquire K 
but K obs , the input vector for a given output target reads 



7>-t retarget 
^in obs 
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We use the setup shown in Figure [T] to record the transmis- 
sion matrix of our system, which is done in approximately 3 
minutes (and require 47V measurements), comparable to the 
time needed to perform an iterative focusing as described in 
fT6l . A speckle spot is chosen as the intended target, and we 
use formula [6] to perform phase conjugation. The results ob- 
tained are shown in Figure [2] 




FIG. 2: Experimental results of focusing, a. initial aspect of the 
output speckle, b. we measure the TM for 256 controlled segments 
and use it to perform phase conjugation, c. norm of the focusing 
operator 0{^ m . d. example of focusing on several points, (the 
insets show intensity profiles along one direction) 

In the single spot focusing experiment (Figure |2]b), we es- 
timate the experimental SNR defined by the ratio of the inten- 
sity at the focus to the mean intensity of the speckle outside. 



The result of the focusing on this target shows a SNR e:E p = 54 
which in par the results of the focusing technique used by 
Vellekoop et al. lfT6l . 

It is important to quantify the effect of both the phase-only 
modulation and the reference speckle used to acquire the TM 
of our complex medium. Indeed, rather then 0^ oc , the effec- 
tive focusing operator includes a normalization to take into 
account the input phase-only modulation. For a single spot 
focusing, 0££ m = K obs Kl orm where k?° rm = ^-/|%|. 
In the representation of the focusing operator 0^. m (Figure 
[2jc), each line corresponds to the expected output intensity 
image resulting while focusing on any point of the CCD's ob- 
servation window. The strong diagonal proves our ability to 
focus light on any of these spots. This is examplified in Fig- 
ure |2]d, which shows focusing on 3 spots. It is worth noting 
that presents weaker sub-diagonals which are due to 

correlations between neighbouring pixels. Finally we would 
like to emphasize that the strong advantage of this focusing 
technique, compared to a sequential algorithm, is that it is not 
necessary to acquire anymore information to change the fo- 
cusing spot or to focus on multiple spots in one step. 

From a theoretical point of view, monochromatic phase 
conjugation focusing through a multiple scattering medium 
has been studied by Derode et al. l22l . The general formula- 



tion for the efficiency of this focusing is SNR ~ N g . 



the 



total number of "information grains", or degrees of freedom. 
In our case, this number is the number of independent pix- 
els (or channels) which is the number N of pixels controlled. 
In our experiment, the use of a reference speckle limits the 
number of degrees of freedom that are controlled. Denoting 
7 the fraction of the SLM used to acquire the TM, and using 
the theoretical formalism developed in [22 ] the effective SNR 
including the effect of the reference writes 

SNR re/ = tt/4 ( 7 7V + (1 - 7)/7) ~ O.877V V7V > 1 

(7) 

Where (1 —7) is the effect of the reference speckle which 
contributes as noise at the focal spot, and the prefactor is due 
to the phase-only modulation l fl"6lL in agreement with a similar 
experiment in acoustics (23). The experimental SNR of the 
focusing shown in Figure [2] b. is 41% of SNR re /, with 7 = 
65%. 

We have shown that the TM of the medium is a powerful 
tool that can be used to focus through a medium, but we want 
to emphasize here that that it extends well beyond focusing. 
An example is the detection of an amplitude or phase object 
placed at the entrance of the scattering sample. Due to the 
spatial reciprocity of wave propagation inside the scattering 
medium, the detection issue is reciprocal to the focusing one. 
Once the matrix is measured, one can retrieve the initial input 
Eobj signal from the output speckle E out = K obs .E ob j with 
the focusing operator 0^ oc . Indeed, the reconstructed image 
reads 

Eimg * Klbs E out = O foc ^E obj = O foc E obj (8) 
Which is almost equivalent to the result obtained in a fo- 



4 



cusinR scheme since O foc as the same statistical properties 
than O^orm ^ ut w i mou t the phase-only effects. In order to 
demonstrate the validity of this approach we performed the 
following detection experiment: instead of placing an object 
between the SLM and the medium, the object is generated by 
the SLM itself using two phase masks E^- and E^-. The 

mask is a plane phase and E^- is obtained by flipping 
the phase of N Q bj pixels from to n. The phase and ampli- 
tude E^ t and E^ t of the output speckles are measured. The 
virtual object E^- — E^. is an amplitude object, recovered 

by calculating E img = K obs (E^ t - E^ 2 J t ). It gives the im- 
ages that would be obtained if an opaque screen drilled with 
Nobj holes of a speckle grain size was placed at the pupil of 
the input microscope objective. The results of this operation 
are shown in Figure [3] for N ob j = 1 and N ob j = 2. It is clear 
that a simple object placed at the input of the sample can be 
reconstructed after propagation, regardless of the randomness 
of the medium. Considering that this imaging method is re- 
ciprocal of the focusing one, the detection is achieved with the 
same SNR. 




a. b. 



FIG. 3: Experimental results of detection of an object placed before 
the scattering sample. Image a. (resp b.) shows the result of one spot 
(two spots) detection for 256 controlled segments on the SLM. 

These experimental focusing and detection experiments 
through a multiple scattering sample are clear demonstration 
that the measured transmission matrix is physical, i.e effec- 
tively links the input optic field to the output one, both in am- 
plitude and phase. Theoretically, this matrix is related to the 
scattering matrix and should be an extremely powerful tool 
for the study of random and complex media. As a first step in 
this direction, we performed the simplest analysis of the TM: 
a singular value decomposition. Indeed, provided that there is 
no ballistic wave in the field at the output of the sample, Ran- 
dom Matrix Theory (RMT) predicts that for those matrices the 
statistical distribution p(X) follows the so-called "quarter cir- 
cle law" where A are the singular values (SV) normalized by 
the total intensity l24l . 

The effect of the diagonal matrix S re f on K is to multiply 
each line by a value Sjj . This changes the variance of the TM 
which affects the SV distribution. To remove this effect, we 
estimate the absolute values of S re f matrix by averaging the 
K obs matrix along every lines. We obtain the diagonal matrix 
S abs defined by sf 3 s = (\k$°\)i ex \ Sjj \ for (\k tJ \) = 1. S abs 
has has no zero on the diagonal thanks to the finite pixel size. 




FIG. 4: Singular values distribution of the experimental transmis- 
sion matrices obtained by averaging over 16 realizations of disorder. 
The solid line is the quarter-circle law predicted for random matri- 
ces. With ■ the matrix filtered to remove the reference amplitude 
contribution and with o the matrix obtained by filtering and remov- 
ing neighboring elements to eliminate inter-element correlations. 

We filter the K obs matrix and obtain : 

Kfu = K x a K x S (9) 

^abs 

Now, is a unitary matrix since it is a diagonal matrix 
of complex number of norm 1, corresponding to the unknown 
phase of the reference speckle. Thus the singular values of 
Kfu are the same as those for K. We finally remove one 
element out of two to remove inter-element correlations be- 
tween nearby pixels. Results of SV distribution on the ob- 
served matrix after filtering and under- sampling are shown in 
Figure [4] the final matrix obtained follows correctly the ex- 
pected quarter-circle law. This distribution differs from the 
bimodal one predicted by lf25l since there is no total energy 
conservation conditions as the transmission matrix measured 
is a small part of the whole scattering matrix 

To conclude, we have demonstrated how the optical trans- 
mission matrix of a multiple scattering media can be mea- 
sured both in phase and amplitude using a spatial light modu- 
lator and a full-field interferometric measurement. We proved 
the validity of the measured TM through focusing and detec- 
tion experiments. Finally, the singular values distribution was 
shown to follow a quarter circle law, characteristic of random 
matrices. We believe that this technique will have many ap- 
plications ranging from the control of wavefield on detection 
and imaging in random and scattering media to very funda- 
mental problems associated with the propagation of wave in 
disordered systems. 
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